Superfluid-Superfluid Phase Transitions in Two-Component Bose System 
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Depending on the Hamiltonian parameters, two-component bosons in an optical lattice can form 
at least three different superfluid phases in which both components participate in the superfiow: a 
(strongly interacting) mixture of two miscible superfluids (2SF), a paired superfluid vacuum (PSF), 
and (at a commensurate total filling factor) the super-counter-fluid state (SCF). We study universal 
properties of the 2SF-PSF and 2SF-SCF quantum phase transitions and show that (i) they can be 
mapped onto each other, and (ii) their universality class is identical to the (d+1 )-dimensional normal- 
superfiuid transition in a single-component liquid. Finite-temperature 2SF-PSF(SCF) transitions 
and the topological properties of 2SF-PSF(SCF) interfaces are also discussed. 
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Recent remarkable success in experimental study of 
ultra-cold atoms in a 3D optical lattice (OL) pj signals 
a major breakthrough in the field of strongly-correlated 
quantum lattice systems. A simple theoretical framework 
which adequately describes physics of atomic gases in OL 
is given by the on-site Bose-Hubbard model U • Its sem- 
inal prediction Q — the superfluid (SF)-Mott-insulator 
(MI) transition — has been unambiguously confirmed 

Realistic experimental perspectives of trapping several 
atomic species in ultra-quantum regime have inspired 
theoretical studies of multi-component systems in OL 
[J HI El I n particular, it has been argued that the 
two-component commensurate mixture of inconvertible 
atoms can be in the super- counter- fluid state (SCF) 
In this state, the net atomic superfluid current is impossi- 
ble, and yet the equal-current flows of two components in 
opposite directions are superfluid. Another intriguing su- 
perfluid groundstate which exists in the two-component 
(with equal particle numbers) Bose system in OL is the 
paired superfluid vacuum (PSF) 0, |3 . Qualitatively, 
this state is equivalent to the superfluid state of two- 
atomic molecules and a BCS superconductor. An impor- 
tant quantitative difference with the BCS theory is that 
bosonic superfluidity exists without pairing correlations 
too, and PSF is always associated with finite intra-species 
interaction. At the moment, it is not clear whether PSF 
can be realized in atomic gases without OL. At the 2SF- 
PSF transition point the pairing interaction is necessarily 
strong, i.e. the scattering length for atoms ready to form 
a pair is of the order of (or larger) than the interatomic 
separation. Under these conditions, metastable atomic 
gases are likely to become unstable from the experimen- 
tal point of view because of very large inelastic cross- 
sections leading to formation of fast tight molecules (not 
to be confused with loose pairs we are discussing here) 
and fast atoms. In OL, this recombination channel is not 
an issue since now the regime of strong/weak interaction 
is controlled by the ratio of the tunnelling constant to the 



strength of the effective on-site interaction, while decay 
rates are still controlled by the one-site physics and are 
not sensitive to tunnelling. 

In this Letter, we discuss universal properties of the 
2SF-PSF and 2SF-SCF phase transitions. First, we prove 
that the two transitions are equivalent to each other by 
establishing mapping between the PSF and SCF phases. 
According to mapping, SCF can be viewed as a "molecu- 
lar" superfluid, where "molecules" consist of particles of 
one component and holes of another component. Corre- 
spondingly, the SCF transition is equivalent to binding 
of two atomic superfluids into PSF. Our main focus is on 
the quantum phase transition. We present strong argu- 
ments that this transition is in the (d + l)-dimensional 
U (1) universality class, and propose an effective (d + 1)- 
dimensional classical model describing it. It allows us 
to relate correlation functions in terms of the original 
bosonic fields to correlators of the U(l) order param- 
eter. In the vicinity of the quantum phase transition 
point, our considerations are naturally generalized to the 
finite-temperature case, predicting the same U(l) uni- 
versality class (but now in d dimensions) for the 2SF- 
PSF(SCF) transition at T > 0. We verify our predictions 
numerically by performing Monte Carlo simulations of a 
3D two-component closed-loop current model of Ref. 
which long-range critical behavior is identical to that of 
a two-component 2D quantum system. Finally, we note 
that the 2SF-PSF(SCF) phase transition preserves the 
"molecular" part of the order parameter and discuss con- 
sequences of this fact for the structure of vortices in hy- 
brid systems containing 2SF-PSF(SCF) interfaces. 

The qualitative equivalence of PSF and SCF phases 
can be understood on the basis of Feynman's represen- 
tation of quantum statistics in terms of particle paths 
(worldlines) in imaginary time t £ [0,(3 — ► oo), where 
[3 is the inverse temperature. In this representation, 
the superfluid groundstate is characterized by worldlines 
forming macroscopic cycles (for brevity, we call them M- 



2 



cycles), when the end of one worldline at r — (3 is the 
beginning of another worldline at time r = 0, and so 
on (partition function worldlines in imaginary time are 
/3-periodic) . The qualitative difference between the 2SF 
and PSF groundstates is that in PSF each A-component 
worldline is bound to some B-componcnt worldline, and 
in the long-range limit there are no free single-component 
worldlines forming independent M-cycles. All M-cycles 
are formed only by pairs of lines, and we arrive at the pic- 
ture of PSF, or molecular superfluid. Less obvious is the 
fact that SCF has the same worldline structure as PSF. 
The key observation is that for integer total filing factor 
one may use a hole representation for one of the com- 
ponents, say, component B. We readily understand that 
the only worldline structure consistent with the absence 
of net superfluid response is when each B-hole world- 
line is bound with some A-particle worldline — this is the 
only possibility of forming M-cycles out of particle-hole 
pairs without having independent single-component M- 
cycles. Macroscopically, bound particle-hole pairs behave 
like new "molecules" with zero particle number charge, 
and their flow is equivalent to the counter-flow of partic- 
ipating components. 

In view of the SCF-PSF equivalence, in what follows 
we discuss PSF only, implying that all results remain 
valid for SCF as well. 

The worldline language presents also a "graphic pic- 
ture" of critical fluctuations driving the PSF-2SF transi- 
tion. Suppose that initially we are deeply in the PSF 
phase. Then each A-line is closely followed by some 
B-line. As the coupling between components becomes 
weaker, bound lines demonstrate local unbinding fluctu- 
ations, see Fig. [Ifa). These fluctuations can be viewed as 
single-colored oriented loops, one part of the loop repre- 
senting, say, an A-line, and the other part representing a 
B-line with the reversed direction, see Fig.QIb). Close to 
the critical point, unbinding loops grow larger and start 
reconnecting with each other (become dense). We assume 
that only large-scale loops are essential for characteriz- 
ing the universality class of the transition; the details 
of short-range behavior are simply determining param- 
eters of the critical action for these loops. The phase 
transition in a system of oriented loops in (d + 1) di- 
mensions (leading to the appearance of macroscopic-size 
loops) is known to describe the SF-MI transition in a 
commensurate system of bosons on a d-dimensional lat- 
tice (see, e.g. 0). In its turn, this transition is equiv- 
alent to the finite-temperature phase transition between 
normal and superfluid states in (d + 1) dimensions 0. 
Hence, the above reasoning suggests a mapping between 
the PSF-2SF and MI-SF transitions, and establishes that 
the PSF-2SF transition is in the universality class of clas- 
sical (d + I)-dimensional U(l) models. 

We now argue that in the long-wave limit our system 
can be also mapped onto a (d + l)-dimensional model 
of two-color classical rotators. This mapping is used to 
have one more argument in favor of the U(l) universality 
class, and to establish important relations between basic 
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FIG. 1: (a) Unbinding fluctuations of two coupled worldlines 
of the two-component system, (b) The single-loop effective 
representation. Worldlines of different species are depicted as 
solid and dotted lines, respectively 

correlation functions; it also provides a natural general- 
ization of our considerations to finite temperatures. 

The presence of lattice commensurability is not cru- 
cial for the PSF-2SF transition since both phases in- 
volved are superfluid. However, it proves convenient to 
formally assume that we are dealing with the double- 
commensurate system — filling factors of both species are 
integer. According to [3j, commensurate d-dimensional 
lattice bosons map in the long-wave limit to a (d + 1)- 
dimensional array of rotators with the Hamiltonian 

ff = - 7 -*,•)> (1) 

<ij> 

where $j € [0, 27r) is the angle of the j-th rotator (or 
phase of the bosonic order-parameter field *£(j) ~ e l * J ) 
and < ■ • • > stands for the nearest neighbor sites on 
a square lattice. In our case, we need three quantum 
fields: ipA and ips for the two components, and the field 
ipp for the pairs. This suggests terms similar to Eq. 
for each of the corresponding three phases. However, 
one also has to account for the terms in the effective 
Hamiltonian converting a bound pair into two atoms and 
vice versa. In terms of the rotator model, this leads 
to a local term oc J2j cos($^ P ' ) — <f>^ 4 ' ) — ^j^) reduc- 
ing the symmetry of the three-color rotor system down 

to U(l) x U(l). This term introduces some (loose) con- 

( p) 

straint on the difference between the phase $^ ; and the 

sum, + . Replacing it with the rigid constraint 

<5>j P ' ) = + we reduce the number of indepen- 

dent variables from three to two — as one could expect 
from the very beginning given the original U(l) x U(X) 
symmetry of our system. As a result, we arrive at the fol- 
lowing Hamiltonian (for simplicity, we assume exchang- 
ing symmetry between the components): 

H = -J2 (7i cos + 72 cos $f + 72 cos ${f>), (2) 

<ij> 
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$ . = *JA) + < , (3) 

where = - and ^f' B) = $f ,B) - <f>f' B) . 
Apart from the 2SF-PSF transition, this model can be 
also used to describe other phase transitions in the doubly 
commensurate system, but not otherwise. 
It is convenient to introduce the phase difference 

<Pi = (*f 5 ~ *$* /2 , (4) 
and to rewrite the Hamiltonian J3J) as (tf ij — tf i — (pj) 
H = - [71 cos $y + 272 cos($y/2) cos . (5) 

[The fields $ and </? describe charge and pseudo-spin de- 
grees of freedom, respectively.] Though the new vari- 
ables, 3>j and tpj, cannot be interpreted as angles of new 
rotators — the configurational space of the original rota- 
tors and is exhausted with, say, $j S [0, 2w) 
and G (— 7r, 7r], while the Hamiltonian (JSj is not 27r- 
periodic with respect to $j — for our purposes it is suf- 
ficient that just ipj can be viewed as a rotator angle. 
Indeed, in both PSF and 2SF, the pair phase variable $j 
is ordered and its local fluctuations are not relevant to 
the criticality of the 2SF-PSF transition. Therefore, we 
may simply set $j = in Eq. J5J which brings us to the 
effective one-component rotor model for tp: 

#2SF-PSF = -272 ^ cosl Pij • ( 6 ) 
<tj> 

The transition thus is the superfluid-normal-fluid tran- 
sition in the (^-channel (which means localization in the 
pseudo-spin sector); the corresponding complex order pa- 
rameter is ip(X.) oc e lv ^ x ', where X is the space-time 
radius- vector treated as a continuous variable in the long- 
wave limit. Given this order parameter and Eq. re- 
lating tp to the original phases (<E>( A ) = — = ^>) 5 we 
immediately find the critical behavior of various correla- 
tion functions 

(V>l(X)v> A (o)> ~ (4(o)V B (X)> ~ {^ A (o)Mx)) 

~(^W(0)). (7) 

Now we note that (d + l)-dimensional model © with 
large but finite size in the r-direction describes the initial 
part of the finite-temperature 2SF-PSF line in the vicin- 
ity of the quantum critical point. We thus establish the 
universality class — U(l) in d dimensions — for the finite- 
temperature second-order 2SF-PSF transition. Since the 
order parameter for the transition is ~ e lip ("molecular" 
order parameter is not critical), we arrive to a rather 
counter-intuitive conclusion that with increasing temper- 
ature the transition is from 2SF to PSF. [Clearly, rela- 
tions take place on the finite-temperature critical line 
as well.] 

This finite-temperature 2SF-PSF transition survives 
even when the two components have slightly different 
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FIG. 2: Different geometries of PSF-2SF interfaces. (See dis- 
cussion in the text.) 

densities and the groundstate is inevitably 2SF (both 
(ipA) and (iPb) are non-zero). Away from the quantum 
critical point, the 2SF-PSF transition can be viewed as 
the superfluid to normal fluid transition in the (dilute) 
sub-system of excessive particles. 

Equation (JSJ) is also useful for understanding the struc- 
ture of vortices across the interface between the PSF and 
2SF phases. Experimentally, interfaces naturally arise in 
inhomogeneous systems (in a confining potential, parti- 
cle densities drops to zero at the boundary, and, e.g. the 
SCF phase, which requires commensurability, may not 
survive at the edge). Suppose one creates a vortex in a 
PSF phase and then adiabatically removes OL and the 
trapping potential to observe the system by the standard 
technique of absorption imaging lla . When the lattice 
potential is turned off, the system will behave as two 
weakly interacting gases. The question now is: Do vor- 
tices in the PSF phase transform (and how) into vortices 
in the resulting weakly interacting system [10|? 

System inhomogeneity implies that at intermediate 
stages of the potential turning off, there will be an in- 
terface similar to one shown in Fig. [2^a) . Ultimately, 
the interface shrinks and disappears with the PSF phase, 
while the topological structure of the 2SF state remains 
the same as it was when the interface existed. To under- 
stand this structure we resort to the rotator model JSJ. 
In both PSF and 2SF the phase field $ is ordered and 
thus the circulation of V<f> does not change across the 
interface. Since the phase field tp is disordered inside the 
PSF phase, there are no topological constraints associ- 
ated with it. Taking into account Q, we arrive at the 
following rule for topological charges: 

J{A) + j(B) =Ij (8) 

where l( A ' B "> and I are integer topological charges in 2SF 
and PSF correspondingly 

t (a,b) = J_ / V $(AB) dl) / = J_Iv$dl. (9) 

2?r Jc 2 27r J Cl 

We see that vortices in PSF always induce vortices in the 
2SF phase fields <&( A > or/and $( B ) , and thus makes it pos- 
sible to observe the circumstantial evidence of the PSF 
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vortex in the absorption image of the weakly interacting 
2SF state. However, the values of topological charges are 
not unambiguously defined. For example, if I = 1 then 
(1^ = = 0) and (1^ = 0,/( B) = 1), are consis- 

tent with Eq. ©, as well as, say, (I (A ^ = 2,I (B '> = -1). 
This implies that particular values of l( A > and 1^ will 
depend on details of the experimental setup determin- 
ing the lowest energy configuration. For example, if the 
two components have different superfluid stiffnesses and, 
initially, there was one vortex in the PSF, then, after cre- 
ating the interface and removing the PSF, the vortex will 
reside in the component with lower stiffness. 

Another interesting geometry is shown in Fig. ffib). 
Using arguments identical to those presented above, we 
see that the only integer topological charge on contour 
C is /. While the sum of integrals for and 1^ still 
equals /, separately they are ill defined on C, because the 
phase if experiences large zero-point fluctuations in the 
PSF region. Suppose, then, that initially there were no 
PSF phase at all, and the topological charges of compo- 
nents were, say, I^ A > = 1 and = 0. Imposing OL to 
create PSF will eliminate quantization for the individual 
phases & A ' B \ while preserving the sum + = 1. 
Accordingly — since no memory about the initial values 
j( b )^j( a ) i s retained — further removal of the OL may 
result in the final 2SF state with 1^ = 0,1^ = 1. 
Similarly to the previous setup, if the two components 
have different superfluid stiffnesses, then, after the cycle 
of switching on and off OL, the circulation will reside in 
the component with lower stiffness. 

There are several options to verify above considera- 
tions numerically. One is a direct simulation of the two- 
component (i-dimensional Bose-Hubbard Hamiltonian at 
very low temperature. However, the universality class 
of the phase transition and the relevant long-wave de- 
scription of critical fluctuations, may be also obtained 
for the (d + l)-dimensional classical lattice model built 
on particle trajectories in discrete imaginary time. One of 
the quantum-to-classical mappings for the Bose-Hubbard 
Hamiltonian — the J-current model — was developed in 



Ref. [9j, and we straightforwardly generalized it to the 
two-component case. Our choice to simulate the classi- 
cal action was motivated only by reasons of numeric effi- 
ciency. Recently developed quantum and classical Worm 
Algorithms do not suffer from critical slowing down , 
but the classical one is superior because of its simplic- 
ity (it was already successfully applied to the disordered 
one-component J-current model 12]). We defer details 
of simulations performed for the d+l = 3 case to a longer 
paper |13j and simply mention here results. The corre- 
lation radius and the correlation function exponents for 
the 2SF-PSF transition agreed with the known values for 
the 3D U(l)-universality class 01 within 1 2% accu- 
racy. We have directly verified that non-trivial relations 
between the correlation functions given by Eq. (|7| hold 
true at the critical point, and deviations from Eq. (JZJ) 
are barely visible even at distances as small as 5 lattice 
constants. We have observed the qualitative prediction 
of the model © about the transition from 2SF to PSF 
with increasing temperature. 

Summarizing, we have shown that two strongly- 
correlated phases of the two-component bosonic 
system — the superfluid state of pairs and the counter- 
flow superfluid — are macroscopically equivalent. We 
have presented arguments supported by results of 
numeric simulations, that quantum phase transitions 
leading to formation of these phases from the state of 
two miscible superfluids are in the universality class 
of superfluid-Mott-insulator transition in a single- 
component bosonic system. The finite-temperature 
2SF-PSF(SCF) transitions are in the universality class 
of a single-component superfluid-normal-fluid transi- 
tion. The proposed two-color rotator model correctly 
describes the critical behavior of various correlators 
and — in the spatially inhomogeneous case — yields a 
simple rule for "sewing" topological defects across the 
boundary between the phases. On the basis of this rule it 
is possible to observe the evidence of the 2SF-PSF(SCF) 
phase transitions even without directly detecting it. 
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